Statistics of static avalanches in a random pinning landscape 
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We study the minimum-energy configuration of a d-dimensional elastic interface in a random 
potential tied to a harmonic spring. As a function of the spring position, the center of mass of the 
interface changes in discrete jumps, also called shocks or "static avalanches" . We obtain analytically 
the distribution of avalanche sizes and its cumulants within an e = 4 — d expansion from a tree 
and 1-loop resummation, using functional renormalization. This is compared with exact numerical 
minimizations of interface energies for random field disorder in d — 2, 3. Connections to the Burgers 
equation and to dynamic avalanches are discussed. 



In numerous systems, the equilibrium or non- 
equilibrium response to perturbations is not smooth and 
involves jumps, avalanches or bursts. In systems at the 
brink of instability, with many metastable states, it is of- 
ten self-organized and critical with power-law tails for the 
probability of large events. This is observed ubiquitously 
in systems with heterogeneities, such as Barkhauscn noise 
and hysteresis in magnets, field response of supercon- 
ductors, contact line of fluids, cracks, granular matter, 
dry friction and earthquakes. Sandpile automata [ij and 
elastic media pinned by quenched disorder [2| have been 
studied as simple models for these phenomena. Rela- 
tions between sandpiles and periodic interface depinning 
and between sandpile models and loop-erased walks 
0] have proved fruitful, especially in d = 2, where 
conformal field theory can be used Q. Despite much ef- 
fort it has proven difficult to obtain analytical results, 
e.g. for the distribution of the size s of avalanches (de- 
fined below), except in mean-field models for sandpiles 
[Toj ] and for random field Ising magnets Q as well as for 
a toy model for avalanches at depinning {2}, which all 
yield P(s) ~ s~ T with r = 3/2. Scaling arguments for 
sandpiles 0, 0] and for depinning fill w ere developed to- 
gether with numerical analysis [3j,ll2j,[l3(. The Functional 
Renormalization Group (FRG) theory for pinned systems 
has led to detailed predictions for e.g. the roughness of 
interfaces but, until now, has failed to describe discon- 
tinuous jump processes 14, 15, lq |. Hence it remains an 
outstanding issue to find a limit where mean field theory 
is valid, prove this, and to develop a controlled field- 
theoretic expansion around it. It should allow to clarify 
the differences between equilibrium and non-equilibrium 
avalanches, recently questioned in a model for magnetic 
hysteresis [17j . 

The aim of this Letter is to provide a first analyti- 
cal calculation of the distribution of avalanche sizes in a 
static, equilibrium setting, using FRG, and to compare 
with numerical calculations. It opens the way to a closely 
related calculation for depinning 18]. As demonstrated 
in our previous work [l^ ] , a model which allows a precise 
FRG treatment and comparison with numerics, both in 



statics and dynamics, consists of an elastic interface in 
a random potential, parameterized by a (scalar) height 
field u[x), and submitted to an external parabolic well, 
i.e. a spring, centered at u = w, 

H[u;w] = J A|[Vii(i)f+y(a:,w(a;))+y[ii(i)-w] 2 . 

(1) 

We are interested in energy minimization as w is varied in 
a given realization of the random potential V(x, u). We 
denote V(w) — mins u r x \\ TL[u; w] the optimal energy and 
u(x; w) the optimal interface position. The force per unit 
volume exerted by the spring is V'(w) = m 2 [w — u(w)], 
where u(w) := L~ d f d d xu(x;w) is the center-of-mass 
position and L d the volume of the system. We study 
the three important universality classes for disorder, with 
short range (random bond, RB), long range (random 
field, RF) or periodic disorder correlator (random pe- 
riodic, RP). Their definitions can be found in 11911 and in 
standard papers on pinning [2(3| and FRG 14 . Il5l . [l6| . 

Although we often use the language of dynamics, one 
should emphasize the difference between the static prob- 
lem studied here, where the interface finds the global en- 
ergy minimum for each w, and the dynamic one, where 
w(t) grows very slowly, and the interface visits a deter- 
ministic sequence of metastable states. In the scaling 
limit to — > 0, on which we focus here, the first case is 
about interface configurations of zero-temperature equi- 
librium, studied in [19| . whereas the latter one is about 
critical depinning, studied in 24, 2^]. In the statics u(w) 
is a (single-valued) function, while it shows some history 
dependence at depinning. Despite these differences, de- 
pinning and statics are close cousins and some differences 
within the FRG are found o nly beyond one loop flit ]. 

As shown previously 0, EE Hfl, [jj, [jj], the optim al 
interface is statistically self-affine with (u(x) — u(0)) 2 ~ 
|x| 2 ^ and a roughness exponent £ which depends on the 
class of disorder, and with a e = 4 — d expansion [lij ]: 
C = e/3 for RF, C = for RP, and C = 0.2083e+0.00686e 2 
for RB (and £ = 2/3ind=l). This holds for scales 
L c < L < L m , where L c is the Larkin length (here of the 
order of the microscopic cutoff) and L m ~ 1/to, the large 



scale cutoff induced by the harmonic well. It is useful to 
picture the interface as a collection of (L/L m ) d regions 
pinned almost independently. 

We found that u(x; w) is an increasing function of w 
which can be decomposed into smooth parts, negligible 
in the scaling limit m — > and jumps (alias "shocks" or 
"static avalanches") as u x (w) — J2iSfO(w — Wi), where 
Si = J d d x Sf is the size of the shock or avalanche la- 
belled i, and 6{x) the unit-step function. The avalanche- 
size distribution, defined from an average over samples, 



J2 S(S - Si)S(w - Wi) = p(S) = Po P(S) 



(2) 



can equivalently be defined from a translational average 
in a given sample. Here P(S) is the normalized size dis- 
tribution and padw the average number of shocks in an 
interval dw. The scaling ansatz 



p(S) = L d m p S- T p{Sm d+ t) 



(3) 



is shown below to hold within the e expansion and verified 
by our numerics. The constraint u'(v) = 1 relates the 
shock rate to the first moment L~ d po(S) = 1, where 
(5™) := J dS S n P(S) denotes the normalized moments. 
It implies for r < 2 the exponent relation p = (2 — r)(d + 
(,) . The distribution is qualitatively different for (i) r < 1 
when a unique scale S m ~ m~( d+ ^ exists, i.e. p(S) = 
S^-piS/Sm), and (ii) 1 < r < 2, where 



P(S) = C T Sv\S/S )- T f(S/S T , 



(4) 



and typical avalanches are of the order of the microscopic 
(UV) cutoff So, while moments (S p ) with p > r — 1 are 
controlled by rare avalanches of size ~ S m , the large-scale 
cutoff. 

As discussed recently [26], shocks for manifolds are 
a natural generalization of shocks in decaying Burg- 
ers turbulence, seen as their d — limit (with, for 
the RF case, t = 1/2 in that limit). The force 
field V'(w) = m 2 [w — u(w)] generalizes the Burg- 
ers velocity field. Its n-point connected correlation 

V'{w 1 )...V'{w n ) = L-^ d (-l) n C ( - n Hw 1 ,...,w n ) 
obeys FRG equations, which generalize the hierarchies 
of Burgers multi-point correlations. Shock-size moments 
can be extracted from their non-analytic part: E.g. the 
second moment is contained in the cusp of the dis- 
order correlator A(w) = C^ 2 \w,0) of the FRG, as 
seen from the relation —m~ 4L L~ d A"(w) — m'(w)m'(0) = 
L~ 2d J2i Sifi( w ~ Wi)6(w) + smooth function, which upon 
integration gives -2A'(0+) = mHS 2 )/{S). This 
generalizes to higher cumulants 26] K^ n \w) := 
m 2n L( n -V d [u(w) - to - u(0)] n ° (the linear cusp of 
measured in [19(, generalizes Kolmogorov's law). Hence 
we can compute the avalanche-size distribution using the 
generating function 



L- d ( e * Ld W«) -«-«(o)] _ i) = Z(X)w + 0{w 2 ) 



I.: 



FIG. 1: Example of a diagram at MF level, as generated by 
Eq. J6) at a = 0. It contains a correction to disorder i.e. 
A'(0 + ) at 1 loop (shaded in gray). 
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for w > 0. It assumes a linear cusp (i.e. a finite density 
of shocks) and generalizes to depinning [3] • 

We have computed the leading non-analyticity of the 
functions and K n (w) 28] from (i) a Legendre trans- 
form of the replicated effective action T computed order 
by order in e; (ii) a direct perturbative expansion without 
replica. The calculation is more involved than usually for 
FRG: the size distribution already at order 0(e ) requires 
a summation of all tree diagrams. The latter could be 
termed mean field, but with the proviso that the scale of 
S involves A'(0 + ) computed to 0(e). Here we compute 
to order 0(e), which amounts to sum all trees and single 



loops; for details see [27|. The main result is that Z{\) 



satisfies a remarkable self-consistent equation to one loop 



Z{\) = A + Z{\) 2 + a y(n + \)T- 2 i n Z(\) 



n>3 



(6) 



where Z(X) 



A'(0+ 

In/{eh), In = J k (k 2 + I)"" 

can graphically be written as 



Z(Am- 4 |A'(0+)|) - A, i r< 



a = -eI 2 m- e A"(Q+). It 



The type of resummed diagrams is presented on figure [T] 
Since a = 0(e), to leading order one solves ((6]) setting 
a = 0. This yields Z(X) = ~(1 - \/l - 4A), identical to 
the generating function of the number of rooted binary 
planar trees with n leaves [29j |. and a size distribution, 
with r = 3/2: 



P(S) 



_ (jV c-1/2 c-3/ 



-S/(4S m ) 



(7) 



This is valid for S ^> So, such that [30( the moments 
with p > \ satisfy (S p )/(S) = UpS^ 1 with a p = 
2 2p ~ 2 7r~ 1 / 2 r(p — i), independent of the non-universal 
small-scale cutoff Sq. Hence the rigorous summation 
of tree diagrams in the FRG yields the same P(S) as 
that of a mean-field toy model for dynamic avalanches 
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FIG. 2: Z(X) for RF, d = 2. The mean field and one loop 
analytical curves are given in the text, A being rescaled so as 
to reproduce the numerically measured second moment (given 
by the curvature of Z(X)). While the MF result differs sub- 
stantially from the numerical measurement, the 1-loop curve, 
i.e. (|9]) with a given by (|8]) setting e = 2 and £i = 1/3, coin- 
cides for all negative A, and almost up to the singularity for 
A > 0. Changing e from 2 to 2.1 or 1.9 already results in a 
visible disagreement for A < 0. 



2] and that of mean-field sandpiles [10]. In addition, 
since the FRG is a first principle method, it predicts 
S m = cm _d ~^ where c = (e/2)|A'(0 + )| is obtained 
from the FRG fixed point for the rescaled correlator 
A(u) = (eJ2)m" e+2 ^A(um'^) and depends on the uni- 
versality class [27j]. Since r > 1 the scale (S) — 5 r " 1 S , ^~ r 
remains undetermined and UV cutoff dependent. Eq. 
seen as a convolution equation for P(S), may allow to put 
the physical picture in [2| on a more rigorous footing. 

To next order in e we solve Eq. ^ which includes 
higher branchings with a universal dimensionless rate 



1 



a 



(l-Ct)e 



(8) 



at the fixed point, where ( = (ie+0(e 2 ) and i n = l/2(n- 
l)(n - 2) in d = 4. It yields 



Z{\) 



l r 

2 



1 - Vl -4A 



4V1 - 4A 



log(l - 4A) x 



x (3A + Vl -4A - 1) - 2(2A + Vl - 4A - 1) 
+0(a 2 ) , (9) 

from which one can calculate the universal ratios: 

2ri- 1 



r n : = (S n+1 )(S n - 1 )(S n )- 2 = 



2n - 3 



(10) 



e nr(n-|) + y¥r(»-i) 

~3 (1 - Cl) (2n-3)T( n -§) +0{€h 

for any real n > 3/2, with (i = 1/3 for RF, £i — for 
RP and Ci = 0.283 for RB. Upon inversion of the Laplace 
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FIG. 3: Numerically computed normalized avalanche distri- 
bution P(S), for random-field disorder and d = 2 (f = 2/3), 
multiplied by S T with r = 1.25 (from (|14[0 to emphasize devi- 
ations from the power law P(S) ~ S~ 1,2l ' . Error-bars are 3a- 
errors for P(S) and the size of the box for S. The solid curve 
is a 1-parameter fit to Eqn. (|11|) . with Sm = 3500, r = 1.25, 
a given by © with = 1/3, e — 2 and the corresponding 
values for B, C given in the text. We use the measured value 
for (S) in (lll[) hence there is no additional free parameter. 
The dashed line is a constant (guide to the eye). Inset: blow- 
up of main plot. The best fit to a pure power law would give 
r = 1.23, and to a power law times exponential r = 1.2. 



transform one finds: 

for S > So, with C = -^y/ira, B = 1 - a(l 
A= 1 + |(2— 37e)«, 7e = 0.577216, and exponents: 

r = l + la= 3 -- 1 -(l-Ci)e + 0(e 2 ) (12) 




* = i-5 = i + i(i-d) e 



(13) 



Note that the decay of large avalanches becomes 
stretched (sub) exponential (in d = for RF, (5 = 3). 



RF 


V2 


T3 


T4, 


mean field 


3 


1.67 


1.4 


d = 3, eq. (10) 


2.33 


1.54 


1.34 


d — 3, numerics 


2.25±.05±.2 


1.48±.04±.14 


1.27±.02±.13 


d = 2, eq. ([10]) 


1.66 


1.42 


1.28 


d — 2, numerics 


1.95±.02±.06 


1.38±.02±.06 


1.21±.02±.06 



TABLE I: Universal amplitude ratios with statistical and sys- 
tematic errors (in this order) for numerics; there is a system- 
atic error since the measured ratios decrease with decreasing 
mass. For d — 2, the decrease which we take as systematic 
error was measured from masses m = 0.025, m 2 = 0.00125, 
and m 2 = 0.000625 (whose values are given). For d = 3, the 
corresponding one is measured for the two smallest masses 
m 2 = 0.0025 and m 2 = 0.00125 (with values from the latter). 
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We note that our result for r agrees to 0(e) with the 
conjecture 



t = 2 - 



(14) 



equivalent to p = 2. It was presented for depinning 11 | 
and for the t s =4/3 exponent of the number of topplings 
in a sandpile model in d = 3 Q , which may be compared 
to CDWs. Since assumptions leading to (fT4")l are not 
rigorous, our first-principle calculation confirms this to 
one loop, and leaves open the possibility of higher-loop 
corrections. Our results are straightforwardly extended 
to the case of LR elasticity [27J ■ 

Exact numerical calculation of minimum-energy inter- 
faces has been performed using the discrete models and 
algorithms as described in [l9[ were details can be found. 
The (corrected) distribution is presented on figure [31 for 
the example of RF in d = 2. We measure 



t = 1.25±0.02 (RF, d=2), r = 1.37±0.03 (RF, d = 3) 

(15) 

This is compatible with eq. (fT4|) . Note that the extra 
stretched-exponential term C in (fTTj) (which could not 
be interpreted as summation of a pre-exponential power) 
leads to a bump which can clearly be seen in the nu- 
merics on figure[31 Finally we have measured (see Fig.QJ 
the distribution of the intervals between successive jumps 
(occuring at positions w = w l ) and found it to be very 
close to a pure exponential. 

To conclude, using Functional RG we have performed 
an expansion around the upper critical dimension to ob- 
tain the avalanche or shock distribution in the statics. 
It compares well with the numerics. Preliminary results 
[lit indicate that the above mean-field and 1-loop results 
also hold for depinning (with the corresponding values 
for C) ; 2-loop calculations are in progress to further check 
the conjecture (fT4"]) and quantify the difference between 
static and dynamic avalanches. 
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FIG. 4: Distribution of intervals between jumps, in units of 
the step-size 8w = 0.002. Error-bars are 3<r-errors. (The last 
5 boxes are overlapping.) 
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